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We show that many-body systems with conserved particle number which have the symmetries 
corresponding to a nonsymmorphic space group have low lying excitations for certain integer values 
of the particle number per unit cell. These results may be interpreted as a generalization of some 
band touching theorems to interacting systems. 



There are relatively few results for non- 
interacting systems that also hold in the pres- 
ence of interactions. Kohn's theorem [l[ and 
Luttinger's theorem are two celebrated ex- 
amples. More recent additions to this list in- 
clude various extensions of the Lieb-Schultz- 
Mattis theorem as applied to number conserv- 
ing systems of fermions or bosons 0, 01 ■ In the 
current work, we establish a technique to ex- 
tend band theoretical results to interacting sys- 
tems and use it to derive interacting counter- 
parts to two well known band touching results 
for nonsymmorphic crystals. The first band the- 
ory result asserts that in systems with certain 
nonsymmorphic space groups, there are specific 
wavevectors at which the electronic eigenstates 
must necessarily be degenerate @; the second 
that for all but two of the nonsymmorphic space 
groups, there are certain directions in reciprocal 
space along which any energy band of the crys- 
tal must necessarily touch another band [1, u\ ■ 
Together, these results apply to all the 157 non- 
symmorphic space groups, which account for 
nearly two thirds of the total number of space 
groups. 

One way to generalize the notion of a Bril- 
louin zone to many-particle systems with inter- 
actions is to consider the phase space associ- 
ated with inserting solenoidal fluxes through the 
holes of a periodic system. Two Hamiltonians 
which differ only by a flux insertion of multiples 
of 2tt through the holes of the torus are equiv- 
alent up to gauge transformations. From this 
perspective, the phase space of flux insertions 
has properties similar to the space of wavevec- 



tors for the single electron problem. Since flux 
insertions do not destroy lattice translational 
invariance, at any point in this phase space, 
the energy eigenstates of a periodic Hamiltonian 
may be labeled by their crystal momenta. The 
energies and the eigenstates associated with a 
particular crystal momentum must form con- 
tinuous functions of the phase space variables, 
which may therefore be regarded as bands. We 
show that for a system with a nonsymmorphic 
space group symmetry and a conserved parti- 
cle number and filling fractions (the number of 
particles per unit cell) equal to certain integers, 
these bands necessarily touch [1] . 

The key implication of these many-body band 
touching results is the existence of low lying ex- 
citations in the limit of large system sizes for 
integer fillings and in the presence of nonsym- 
morphic space group symmetries. We can es- 
tablish a lower bound on the number of these 
excitations, which depends on the filling frac- 
tion and the nonsymmorphic space group. 

The arguments that we are about to present 
are quite general and can be extended to a vari- 
ety of systems, some of which will be discussed 
at a later stage. It is simplest, however, to 
first consider a system of spinless fermions or 
bosons in 3d with a conserved particle num- 
ber which is invariant under translations by the 
primitive lattice vectors, ai,a2 and as. We will 
also assume that the filling fraction is some in- 
teger V. We impose periodic boundary con- 
ditions (TaJ^ = (TaJ^ = (ra3)^ = 1 for 
some integer, N , where is the operator that 
corresponds to a lattice translation by r. Let 
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bi,b2,b3 be a basis for the reciprocal lattice 
such that b^.aj = 27r(5y . 

We will assunie that the system is also invari- 
ant under all the symmetry operations corre- 
sponding to some nonsymmorphic space group, 
Gi. 

Consider matrix representations of the group 
of k [6,] for which the matrix corresponding 
to a translation by an arbitrary lattice vector, 
r has the form e^^ '^I where / is the identity 
matrix. Let nk be the dimension of the low- 
est dimensional irreducible representation of the 
group of k with this property. Suppose further 
that there is some ks for which rik^ is strictly 
greater than 1. Our aim is to show that for 
groups which have such a point and for which 
some additional conditions hold, there is a finite 
set of points in the phase space of solenoidal fiux 
insertions at one of which the ground state (or 
any other eigenstate) is necessarily degenerate. 

Owing to a similarity in structure, the space 
of flux insertions is conveniently parametrized 
by vectors in reciprocal space. A solenoidal 
flux insertion of 6i , 02 and 6*3 through the holes 
of the 3 torus corresponding to the directions 
ai , a.2 and a.^ may be called the flux insertions 
associated with the wavevector k = dihi + 
6'2b2 -I- 6'3b3. Let -ff(k) be the Hamiltonian 
which includes the effects of flux insertion as- 
sociated with k. Rather than using periodic 
boundary conditions with H(k), we may work 
with the original Hamiltonian, H and the gener- 
alized periodic boundary conditions associated 
with k which may be defined as: (TaJ^lV') = 

e^«MV),(ra.r = e^''\^),{T^,r = e^'^l^) for 
single particle states 

For any wavevector, k, the boundary condi- 
tions associated with that wavevector are com- 
patible with the group, of that wavevector. 
To see this, let g be some element of G^ and Rg 
a proper or improper rotation associated with g 
such that tn^^g = gtr for any lattice translation, 
ir in the space group, G. Then, if lip) is some 
single particle state which satisfies the general- 
ized boundary conditions associated with k, and 



D(g) the unitary operator corresponding to g, 

^g.(JVa.).tq,^(^)|^^ (1) 

The last equality for i = 1,2,3 implies that 
D{g)\il)) also satisfies the boundary conditions 
associated with k. 

Now, if we define the twist operators, U{k.): 

where Uj. is the fermion number operator, 
and if is an eigenstate of H{k) with energy 
E (with periodic boundary conditions), then 
[/(k)!?/;) is an eigenstate of H with the same 
energy E with the generalized periodic bound- 
ary conditions associated with k. The twist op- 
erator does not commute with translations and 
other space group operations. Instead, we have 
for any Np particle state, I?/;), 

T,,U{\<i)\^) = e*(^''/^)(''-'^')[/(k)rr-|i^) (3) 

and if D{g) is the unitary matrix correspond- 
ing to a space group operation g consisting of a 
(proper or improper) rotation R followed by a 
translation by r and if i?k = k, then 

D[g)U{\i)\il:) =e'(^''/^)(''-")C/(k)Z?(g)|tA) (4) 

Suppose, further, that there is a set of in- 
tegers, n such that n(ks — ko) is a reciprocal 
lattice vector for any ko such that the group, 
G*' of k = ko is the entire space group. Let 
p be the smallest member of this set and let 
k' = ks — ko. Suppose now that the ground 
state |(/)) of the Hamiltonian, H with no fluxes 
has a momentum, ko such that the group of 
ko is the full space group, G. From our previ- 
ous discussion, we know that there must be a 
continuous band of states, one for each k with 
the same momentum ko. Let |(/'(k)) denote the 
state which belongs to this band and is an eigen- 
state of i?(k). 
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For a system where N is chosen to be 
coprnne with p, consider the set of states 
{{/(k')|0(k')), {/(2k')|0(2k')) • • • , Uipk'Mipk'))} 
From Eq. (jS]), we deduce that the set of mo- 
menta of these states consists of the wavevectors 
ko + mvN^'k' for m = 1,2,. ..p. If v is also 
coprime to p, then the set must contain the 
wavevector kg. There must therefore be an 
Q <l <p such that ko + IvN'^k' = k^. Then, 
since n\^^ > 1, this statement is equivalent 
to the statement that there must be some 
Q < I < p for which the Hamiltonian with 
boundary conditions associated with Ik' or 
equivalently, the Hamihonian, H{lk') must 
necessarily have a degeneracy. If the group 
of kg (the momentum of the ground state) is 
not G, then for some space group element, 5, 
D{g)\(j)) must also be an eigenstate of H and 
has a different crystal momentum and is hence 
orthogonal to Thus in this case there is a 
degeneracy of the Hamiltonian Hilk) for / = p. 
These arguments establish many-body band 
touchings in more than a hundred nonsymmor- 
phic space groups including the space group of 
the diamond lattice, and the two exceptional 
nonsymmorphic space groups which do not 
have a nonsymmorphic element. 

We now present another set of arguments 
based on the monodromy of the irreducible rep- 
resentations of nonsymmorphic space groups [6|, 
Q • An element 5 of G is called nonsymmorphic 
if there is no lattice translation tr such that gtr 
leaves some point in space fixed. For any non- 
symmorphic element g of the space group G, 
there is some smallest integer Ug such that (/"f is 
a translation by some lattice vector, Yg. Let kg 
be the smallest reciprocal lattice vector in the 
direction of r^. Let G'g be the subgroup of the 
space group generated by the primitive trans- 
lations and g. Since G'g is a subgroup of the 
group, G*^ of kg, the boundary conditions asso- 
ciated with any vector of the form tkg (where 
t is some real number) are compatible with the 
symmetry operations in G'g . Thus the Hamilto- 
nian, H{tkg) with solenoidal fluxes associated 
with tkg is also invariant under these symme- 



try operations. 

Let T be the set of crystal momenta k such 
that G'g is a subgroup of the group G^ of k. 
Let be some state in the ground state man- 
ifold of the system with no fluxes, with a crys- 
tal momentum ko G T. We can then always 
choose |</>) such that it is associated with a one 
dimensional irreducible representation, of 
G'g. By continuity, for all t, there must then 
be an eigenstate \(j){t)) oi H{tkg) with the same 
momentum kg and associated with the same ir- 
reducible representation, Z3^^ of G'g. If G'g is a 
subgroup of G^ of k, then there are Ug inequiva- 
lent one dimensional irreducible representations 
of G'g associated with k: 

D^{g) = g^(k-rg/ng+p27i-/nj,) j-g-j 

where p is an integer which takes the values 
l,2---ng. 

As a consequence of Eqs. ^ and (|4]), the 
state U {tkg)\<l)[t)) which is an eigenstate of the 
Hamiltonian with the boundary conditions as- 
sociated with tkg has a momentum f{t) = kg -I- 
vN'^tkg and transforms per the irreducible rep- 
resentation, D^f(fy Since kg is a reciprocal lat- 
tice vector, the boundary conditions associated 
with kg are equivalent to those associated with 
0. Further, since /(I) is a wavevector equiva- 
lent to kg, -0^(1) must be the same as Z?^^ for 
some Q < p' < Ug. In fact, p' is easily seen to be 
{p + vN'^ikg.Tg)) mod rig. It has been shown 
that for every element q o f the type considered 
here, e*('<^« '"!')/"3 7^ 1 This implies that 

there are some z/, N for which p' ^ p mod n. 
For these values of v and N , the inequivalence 
of the representations associated with and 
D^f(^^-^ implies that the kets 10) and U {kg)\(j){kg)) 
are orthogonal. This in turn, means that the 
band associated with must necessarily touch 
some other bands in the space of solenoidal flux 
insertions. 

We may repeat this process of flux insertion 
starting at each new step with the state \(t>') 
generated in the previous step to produce a set 
of orthogonal states, whose number is the order 
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of {i/N^Ckg.rg)) mod ng in the group of inte- 
gers modulo Ug. The maximum possible value 
of this number for a given v and k^.r^ is ng/ fig 
where /3g is the greatest common divisor of Ug 
and v{\<ig.Vg). In order to generate this many 
orthogonal states, we need N to be coprime to 
ng/fi. 

For each of the nonsymmorphic elements g of 
G 11| we can use the same technique to gen- 
erate connected orthogonal states. With an N 
which is coprime to the set of all Ug/ j3g, we get 
the maximum possible lower bound on the num- 
ber of orthogonal states that can be obtained in 
this way, which is the lowest common multiple 
of all of the numbers Ug/ (3g. 

If the momentum kg of |(/)) is not in T, let 
q be the smallest integer such that 5^ G G^^' . 
Then the states \(t)'a) = g''\(t>) for < a < q 
are orthogonal and necessarily degenerate. For 
each of these q states, we can, in general, 
generate additional orthogonal states which cor- 
respond to the element g' — g'^ . The number of 
such states is, using the formula derived above, 
Ug/q"/ where 7 is the greatest common divisor 
of Ug/q and i^(kg.rg). This number is bounded 
below by ng/q/3g, thus the total number of con- 
nected orthogonal states associated with g is 
greater than or equal to ng//3 and our previ- 
ous bound on the total number of orthogonal 
states still holds. 

Let us now consider a system of particles with 
spin. We assume the system is invariant under 
a set of unitary transformations corresponding 
to the double group of some space group, G. We 
note that all the representations of the double 
group of G also satisfy Eqs. ([3]) and (H]). When 
i'{2S)N^ is even, where S is the spin of the par- 
ticles, both our previous arguments which es- 
tablish connectivity of bands apply in entirety 
since all eigenstate manifolds can be decom- 
posed into invariant subspaces associated with 
irreducible representations of the single group. 
When iy{2S)N^ is odd, all the eigenstate man- 
ifolds can be decomposed into invariant sub- 
spaces of the additional irreducible representa- 
tions of the double groups. 



We are thus led to consider the additional ir- 
reducible representations of the double group 
of G^ such that the matrix corresponding to a 
translation by any lattice vector r has the form 
e*"^ ""/ with / being the identity matrix. Suppose 
there is some point k^ where the minimum di- 
mension of the additional irreducible represen- 
tations of the above form of the double group 
is greater than 1 1^. We assume that, in addi- 
tion, for each kg such that the group, G*^ of ko 
is the entire space group, there is some small- 
est integer p such that p(k^ — kg) is a recipro- 
cal lattice vector. The remainder of our first 
band connectivity argument goes through un- 
changed on substituting kg by k^. Thus, for 
appropriate values of v and N, we may deduce 
that there is a finite set of points in the space of 
flux insertions at one of which the energy bands 
must necessarily touch. From the tables of irre- 
ducible representations of space groups, we find 
that these arguments apply to a large class of 
space groups (more than a hundred). 

To apply our second argument (when 
v{2S)N^ is odd), we note that the additional 
irreducible representations of the double group 
of G'g which involve a single wavevector ki , such 
that G'g is a subgroup of the double group of ki , 
are of the form : 



KM) 



„j(ki.rg/rig + (2p+l)7r/ng) 



(6) 



for < p < Ug. Noting that these irre- 
ducible representations have the same mon- 
odromy properties as the irreducible representa- 
tions of the single groups, we conclude that the 
number of connected orthogonal states that can 
be generated through flux insertion remains the 
same as in the case of particles without spin, 
for the same values of u and TV for the same 
operation, g. 

If the system is invariant under just the space 
group transformations of some group G without 
the accompanying rotations in spin space [13[, 
it can in some cases be mapped to a problem 
with a different filling fraction. Eigenstates of 
spin 1/2 particles in a system with total Sz and 
particle number conservation can be mapped to 
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a two species bosonic model which gives a fill- 
ing factor of 1/2 — m per species. A quantum 
spin system with total conserved magnetization, 
S'fot can be mapped to an interacting boson sys- 
tem with S — m particles per spin, where S is 
the spin quantum number and m is the average 
magnetization per spin. Then, in the presence 
of space group symmetries, a fully gapped phase 
is not possible even for certain integer values of 
n{S — m) where n is the number of particles 
per unit cell. Any plateaus at the correspond- 
ing values of magnetization, if they exist, must 
be "exceptional" [i,[l3. 

What we have so far shown is that under cer- 
tain conditions, the energy bands of the many- 
body systems with integer filling fractions in the 
phase space of solenoidal fluxes must necessarily 
touch. This however necessarily means that for 
these systems, with some additional constraints 
on the locality of the Hamiltonian, there are low 
lying energy states in the appropriate thermo- 
dynamic limit using the same arguments that 
were used by Oshikawa and Hastings for the 
fractional filling case 0, 0, Essentially, a 
gap to local excitations in a system with in- 
teractions that are sufficiently local ensures an 
insensitivity to boundary conditions. This in- 
sensitivity to boundary conditions can be com- 
bined with our arguments to show that for all 
the nonsymmorphic space groups there is nec- 
essarily at least one (and in many cases, more 
than one) low lying excitation at zero flux whose 
gap from the ground state can be bounded by a 
number that goes to zero in the thermodynamic 
limit [is!] . The vanishing of the matrix elements 
of all local operators between the different low 
lying states would then imply the existence of 
topological order. Failing this, the system must 
have some sort of symmetry breaking. 

While our discussion has been in 3d, the same 
arguments can also be applied (with similar re- 
strictions) to nonsymmorphic crystals in 2d to 
show that for certain integer filling fractions and 
appropriate system sizes, there are necessarily 
low lying excitations. Extensions to higher di- 
mensions is left for future work. 
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